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OO , Abstract 

An updated analysis of all available neutrino oscillation evidence in Solar 
experiments including the latest SNO data is presented. Predictions for total 
rates and day-night asymmetry in Borexino are calculated. Our analysis features 
the use of exhaustive computation of the neutrino oscillation probabilities and 
Q_i' of an improved statistical treatment. 

^ • In the framework of two neutrino oscillations we conclude that the best fit 

D . to the data is obtained in the LMA region with parameters (Am^,tan^0) = 

(5.4 X 10"^ eV^O.38), (xLn/^ = 0-81, n = 38 degrees of freedom). Although 
less favored, solutions in the LOW and VAC regions are still possible with a 
^ ' reasonable statistical significance. The best possible solution in the SMA region 

' gets a maximum statistical significance as low as ~ 3%. 

We study the implications of these results for the prospects of Borexino and 
the possibility of discriminating between the different solutions. The expected 
normalized Borexino signal is 0.63 at the best fit LMA solution, where the DN 
asymmetry is negligible (approximately 10~^). In the LOW region the signal 
is in the range ~ 0.6 — 0.7 at 90% confidence level while the asymmetry is 
~ 1 — 20%. As a consequence, the combined Borexino measurements of the total 
event rate with an error below ±5 — 10% and day-night total rate asymmetry 
with a precision comparable to the one of SuperKamiokande will have a strong 
chance of selecting or at least strongly favoring one of the Solar neutrino solutions 
provided by present data. 
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1 Introduction 



The Solar neutrino problem p]-p!0| has been defined, for example, as the difference 
between the neutrino fiux measured by a variety of Solar neutrino experiments and 
the predictions of the Standard Solar Model (SSM). The explanation of this difference 
by the neutrino oscillation hypothesis has been boosted by the SNO demonstration 
at the 3a level of the appearance at the Earth of Sun-directed muon and tau neutrino 
beams |11| and the first determination of the total neutrino fiux generated by the 
Sun [|T^. The agreement of this fiux with the expectations implies as a by-product the 
confirmation of the validity of the SSM |T^-|15|. 

At present the SNO experiment measures the '^B Solar neutrinos via the reac- 
tions [|16]-[g: 1) Charged Current (CC): Ue + d ^ 2p + e- , 2) Elastic Scattering (ES): 
+ + e~. The first reaction is sensitive exclusively to electron neutrinos. 

The second, the same as the one used in SuperKamiokande (SK), is instead sensitive, 
with different efficiencies, to all flavors. The flrst results by SNO on Solar neutrinos 
(1169 neutrino events collected during the flrst 240 days of data, Ref. |p.2||) conflrm pre- 
vious evidence from SK and other experiments |jl|, |^, ^, |20| . The Solar neutrino flux 
measured via the observed CC reaction rate, (pcc'^ = 1-75 ± 0.07 ± 0.12 x 10^ cm~^ 
s~^, is lower than the two existing ES reaction rate measurements: the more precise 
at SK 0|f = 2.32 ± 0.03 ± 0.08 x 10^ cm'^ s'^ and the one measured by SNO itself 
^SNO ^ 2.39 ± 0.34 ± 0.15 x 10^ cm-^ s'\ Comparison of the and and es- 

pecially (pcc'^ and cj)^^ can be considered as the flrst direct evidence for Solar neutrino 
oscillations which cause the appearance of sun-directed muon and tau neutrinos at the 
Earth. When considering only SNO results, it is advantageous to use the experimental 
ratio of CC to ES measurements because of the cancellation of systematic errors. The 
experimental value is 

R'^^P = (j)^cc^ / (t)%^s^ = 0.748 ± 0.13 

which is 2o" away from the no-oscillation expectation R^^"^ = 1. Recall that the sta- 
tistical signiflcance of the difference (pcc'^ ~ 4'%^^ is clearly lower (~ lAa) |jl2[. If we 
use the ES-SK result, the no-oscillation hypothesis is excluded at ~ 3.3(t ^ 



Moreover the allowed oscillations to only sterile neutrinos is excluded at the same level 
(~ 3o"). Although less favored, the oscillations to both active and sterile neutrinos are 
still allowed 

Although of qualitative importance, the SNO results by themselves, due to their 
large error bars, are still far from constraining the neutrino oscillation parameters in 
a signiflcant way. For this purpose we continue to need the data of all other Solar 
neutrino experiments |l], ^, ^, 13^, |3l]] . The aim of this work is to present an up-to-date 
analysis of all available Solar neutrino evidence including the latest CC-SNO data in 
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a two-neutrino framework. The predictions for day and night signal rates in the near 
future real-time ^Be sensitive Borexino experiment [p^-|35[] are also calculated. Two 
important characteristics of this paper are the use a thorough numerical computation of 
the neutrino oscillation probabilities and the use of an improved statistical analysis. 
The latter is explained in detail in the next section. In summary we have followed the 
SK standard procedure as close as possible allowing for a better treatment of correlated 
uncertainties. 

As we show below, the best fit to the full data set including global rates and SK 
energy spectrum is obtained in the Large Mixing Angle (LMA) region. Although less 
favored, solutions in the Low mass (LOW) and Vacuum (VAC) regions are possible 
with a reasonable statistical significance. A best fit solution in the Small Mixing Angle 
(SMA) region is strongly disfavored. 

This work is organized as follows. In Section |^ we present a summary of the 
methods that are used: details of the computation of neutrino oscillation probabilities 
in the Sun and in the Earth, scattering cross sections and calculation of the signal at 
the detectors. We continue with a detailed account of the statistical methods. In 
Section ^ we present the results of different analyses: first considering total event rates 
only and then including the SK energy spectrum information. In Section ^ we present 
our Borexino predictions for the best fit solutions obtained in the previous sections. 
Finally in Section ^ we summarize and draw conclusions. 



2 Methods 

2.1 Neutrino oscillations and expected signals 

In general, our determination of neutrino oscillations in Solar and Earth matter and 
of the expected signal in each experiment follows the standard methods found in the 
literature Nonetheless, this work differs from previous ones in the treatment of the 
computation of the neutrino transition probabilities. We completely solve numerically 
the neutrino equations of evolution for all the oscillation parameter space. Although 
other alternatives, as the use of semi-analytical expressions in portions of the (Am^, 9) 
plane, could be previously justified in terms of economy of resources, at this stage a full 
simulation is an affordable and better option. This is specially true in the calculation 
of the smearing of the neutrino probability as a function of the neutrino production 
point when introducing non radial propagation. However no suprises are found: our 
numerical results are in overall good agreement with the semi-analytical approach. 

The survival probabilities for an electron neutrino, produced in the Sun, to arrive 
at the Earth are calculated in three steps. The propagation from the production point 
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to Sun's surface is computed numerically in all the parameter range using the electron 
number density given by the BPB2001 model [111] averaging over the production 



point. The propagation in vacuum from the Sun surface to the Earth is computed 
analytically. The averaging over the annual variation of the orbit is also exactly per- 
formed using simple Bessel functions. To take the Earth matter effects into account, we 
adopt a spherical model of the Earth density and chemical composition. In this model, 



the Earth is divided in eleven radial density zones |]3^, in each of which a polynomial 
interpolation is used to obtain the electron density. The glueing of the neutrino prop- 
agation in the three different regions is performed exactly using an evolution operator 
formalism [Q. The final survival probabilities are obtained from the corresponding 
(non-pure) density matrices built from the evolution operators in each of these three 
regions. The night quantities are obtained using appropriate weights which depend on 
the neutrino impact parameter and the sagitta distance from neutrino trajectory to 
the Earth center, for each detector's geographical location. 

The expected signal in each detector is obtained by convoluting neutrino fluxes, os- 
cillation probabilities, neutrino cross sections and detector energy response functions. 
We have used neutrino-electron elastic cross sections which include radiative correc- 



tions |38, 39|. Neutrino cross sections on deuterium needed for the computation of the 



CC-SNO measurements are taken from [^^^ 



Detector effects are summarized by the respective response functions, obtained 
by taking into account both the energy resolution and the detector efficiency. The 



resolution function for CC-SNO is that given in [|T2|. We obtained the energy resolution 
function for SK using the data presented in [^l]-^ . The effective threshold efficiencies, 
which take into account the live time for each experimental period, are incorporated 
into our simulation program. They are obtained from |^ . 



2.2 The calculations: definitions and procedures 

The statistical significance of the neutrino oscillation hypothesis is tested with a stan- 
dard method which we explain in some detail here. 



In the most simple case, the analysis of global rates presented in Section |3.1| , the 
definition of the function is the following: 



cxp\ 



where cr^ is the full covariance matrix made up of two terms, cr^ = cr^^c + o^cor- 



The diagonal matrix 0"^^^, contains the theoretical, statistical and uncorrelated errors 
while cr^oj. contains the correlated systematic uncertainties. The i^**^'''^^ are vectors 
containing the theoretical and experimental data normalized to the SSM expectations. 
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The length of these vectors is 3 or 4 depending on whether the CC-SNO experiment 
is included or not: 



^tKexp ^ gth,e.p I gSSM ^ ^ ^ ( CC-SNO ) . 

The index i denotes the different Solar experiments: Chlorine (CI), Gallium (Ga), 
SuperKamiokande (SK) and Charged Current SNO ( CC-SNO ). 

The correlation matrices, whether including the SNO experiment or not, have been 



computed using standard techniques |44 



To test a particular oscillation hypothesis (Am^,tan^^^) against the parameters 
of the best fit, we perform a minimization of the Xgiob ^ function of the oscilla- 
tion parameters. A point in parameter space (Am^,tan^^) is allowed if the glob- 
ally subtracted Xgiob fulfills the condition XgiQ|3(Am^, — xLin < Xn{CL). Where 
X^=2(90%, 95%, ...) = 4.60, 5.99, ... are the n = 2 degrees of freedom quantiles. 

A more elaborated statistical analysis is necessary once one introduces the SK 
energy spectrum into the analysis (Section |3.2D . 



The procedure assumed by the SK collaboration, in order to present its own analysis 
and obtain the tables of bin correlated errors, is that of using the summation = X] 



where each term is of the form [^, ^ 

X^ - am.,)R'l'~Rr, (2) 

and where the expression 

mor) ~ ^ . (3) 

-L ~r ^z.cor^cor 

is the response function for the correlated error in the z^'^-energy bin. The CTj^cor are bin- 
correlated uncertainties and a is an overall flux normalization factor. The correlation 
parameter 5cor is arbitrary and determined in the minimization process together with 
the rest of oscillation parameters. 



SuperKamiokande calculation of tables for correlated errors (as Table 111 in Ref. 
and previous publications) and the prescription given respectively by Eqs. (|^) and (|^) 
are closely related. In this work we follow this definition of the x^? however for clarity we 
write the expressions in a notation which renders the comparison with other published 
works simple. It can be easily shown that the expression to be presented below leads 
to the SK expressions 0. 

Therefore, for the analysis of the full set of data including SK energy spectrum, 
we consider a function which is sum of two quantities = Xgiob + Xspec- The first 
^We have not included a term logdet cr^((5cor) ^ log(l + (5cor) whose effect is neghgible 
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one is given by Eq. (|1]) considering only C/, Ga and CC-SNO total rates, while for the 
second term we adopt the definition 

A subindex i = d,n corresponding to separated day and night quantities is understood 
for any R vector and a matrix. We have introduced the fiux normalization factor a and 
the correlation parameter 5cor- The complete variance matrix is not a constant quantity. 
It is obtained by combining the statistical variances with systematic uncertainties and 
dependent on this correlation parameter. 

The corresponding terms in Eq. (Q) are 

Xa _9 ' Xs 



The central value a"^ = 1 corresponds to the BPB 2001 '^B fiux. The uncertainty on a 
comes from the ^{4% uncertainty in the ^B fiux and we make an average of the two-side 



errors taking a^^ ± (Jq = 1.00 ± 0.17 |T^. The correlation parameter is assumed to be 
constrained to vary in a gaussian way with an error corresponding to as and we take 
S'^ = 1. 

"cot 

In this work we consider two cases. In the first case the fiux normalization is taken 
to be free (we allow for a very large cTq, ^ 00 in the Xa term). As a second possibility 
we consider that both quantities a and 5cor are constrained. Note that the unwanted 
consequence of eliminating the last term x's would be to obtain running-away solutions 
for the optimal value of this parameter. 

The summation now contains 41 bins in total: 3 from the global rates (all the 
experiments except SK) and 2x 19 bins for the SK day and night spectrums. The full 
correlation matrix is defined by blocks. The 3x3 block corresponding to the global rates 
is defined as above. For each day and night spectrum the corresponding 19x19 block 
correlation matrices are conservatively constructed assuming full correlation among 
energy bins [|. The components of the variance matrix are 

i'^cor) ij ~ '^i.expCj.cxp + 0"i,calC'"jf,cal) 
{'-'uncJii ~ "^i.stat "T "^i 



where cxstat are the statistic errors and the quantities o"oxp, o"cai, func are respectively 
the bin-correlated experimental, the spectrum calculation uncertainties and the bin- 
uncorrelated ones. 



^The introduction of the parameter Scot is equivalent to the relaxation of this condition during 
minimization. 
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Notice that in the above sum we do not simultaneously include the global ratio for 
SK and the partial spectrum bins. We ignore any correlation between the spectrum 
information and the global rates of all experiments except SK itself. Finally, a remark 
is in order. In this case the defined covariance matrix is dependent on one of the fitting 
parameters. One does not, however have to add any correction (i.e. log deter terms) to 
the Ax^ expression as long as we sit at 5™™. To test a particular oscillation hypoth- 
esis and obtain allowed regions in parameter space we perform a minimization of the 
four dimensional function x^(Am^, tan^ 9, a, 5cor)- The minimization of the expression 
with respect to a and 5cor is done analytically. The subsequent minimization in the 
(Am^, tan^ 6) plane is numerical. For 5cor = ^Tmi ^ = «min, a given point in the oscilla- 
tion parameter space is allowed if the globally subtracted quantity fulfills the condition 
Ax' = x'(Am2, 9) - xLn < xl{CL). Where xlU^^%, 95%, ...) = 7.78, 9.4, ... are the 
quantiles for four degrees of freedom. 



3 Results 



We have used data on the total event rates measured at chlorine Homestake experiment, 
at the galhum experiments SAGE GNO H and GALLEX and at the 

water and heavy-water SK (live time 1258 days) Q and SNO {ES,CC 240 days |T! 



experiments (see Table (|l]) for an explicit list of results and references). 

For the purposes of this work it is enough to summarize all the gallium experiments 
in one single quantity by taking the weighted average of their rates. 

In addition, for SK we use the available information for the day and night energy 
spectrums. For each of the day and night cases, this spectrum information contains 18 
total recoil-electron energy bins of width 0.5 MeV in the range 5 to 14 MeV and an 
additional bin spanning the remaining 14-20 MeV range. The data and errors for indi- 



vidual energy bins for SK spectrum has been obtained from Ref. [^. The information 
from other SK results as the global day night asymmetry is already contained to a 
large extent in the previous quantities and does not change the results to be presented 
on continuation [H, |H . 



3.1 Global rate analysis 

The analysis of the global rates of the four experiments SK, Chlorine, Gallium and 
CC-SNO is the simplest possibility but nonetheless it reveals important trends of the 
solutions to the Solar neutrino problem. It illustrates how the CC-SNO data impose 
additional constraints on the parameter space and in particular how the SMA solution 
loses its statistical significance in favor of the LMA and LOW regions, becoming allowed 



7 



only at marginal ^ 3a confidence levels. 

We present results both with and without the inclusion of the CC-SNO global rate. 
In the analysis of the global rates of the SK, CI, Ga ( CC-SNO ) experiments one has 
two free parameters Am? and tan^ 9 and 3 (4) experimental quantities, therefore the 
effective number of degrees of freedom (d.o.f.) is 1 (2). 

On the left-hand-side of Tables we present the best fit parameters or local 

minima obtained from the minimization of the function given in Eq. Also 
shown are the values of Xmin P^r degree of freedom (x'^/n) and the goodness of fit 
(g.o.f.) or significance level of each point (definition of SL as in Ref. [^). On the 
right-hand-side of the same tables we show the deviations (minimization residuals) 
from the expected values for all the total event rates. We also include, for further 
reference, the deviations corresponding to the SK day and night energy spectrums 
and global day-night asymmetry, although these quantities are not incorporated in 
the minimization. We first consider the case in which the CC-SNO global rate 
is ignored. The absolute minimum is located at the SMA region (Am^, tan^ 6')'^^^"^ = 
(7.8 X 10"^ eV^ 1.10 X lO'^), with xLJ^ = 1-06 (d.o.f. n = 3 - 2), the significance 
level for this point g.o.f. is 29%. Closely we find the twin minima situated in the VAC 
region with parameters (Am^ tan^ 9)^^^ = (8.4 x 10"^^ eV^ 3.37(0.27)), which receive 
a significance level only slightly smaller than the SMA mininimum. The next local 
minima are situated in the LMA region (Am^, tan^ 61)^*^^ = (1.23 x 10~^ eV^O.34), 
with Xmin/^ — 1-89 which corresponds to a still reasonable good fit g.o.f. is 17% and in 
the LOW region (Am^, tan^ ^)^*^^ = (1.11 x 10"^ eV^0.75), with x'Ljn = 3.62. For 
comparison, much-less significant results are obtained when sterile neutrino oscillations 
are considered. 

When one introduces the CC-SNO total event rate in the analysis the relative 
order of the local minima changes. The absolute minimum is now located at the 
LMA region: (Am^, tan^ ^)^^^ = (3.12 x 10"^ eV^0.40). The significance level of 
the absolute minimum is clearly worse than in the previous case Xmm/''^ — 1-61 (d.o.f. 
n=4-2) which corresponds to g.o.f. = 21%. 

The significance level of the SMA region which follows is no better than Xmin/''^ ~ 
2.40, g.o.f. is 8.8%, value which is obtained at (Am^, tan^ e)^^^^ = (7.3x 10"^ eV^ 1.40x 
10~^). After that, we have the pair of minima at the VAC region (Am^, tan^ 6')^'^'^ = 
(9.7 X 10^11 eV^ 4.00(0.26)), with Xmin/^ = 3.2 - 3.1. Finally, we find in the table a 
minimum in the LOW closely situated to the one of the previous case. 

In Figs. (P we present graphically our results respectively before and after including 
the global CC-SNO rate. In the plots one can see the regions which are allowed at 
90, 95, 99 and 99.743% confidence levels. The main difference between the plots is an 
overall reduction of the extent of the allowed area at any CL. The large contiguous areas 
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at 3 and 4a which are present along tan^ ^ ~ 1 in the first case are converted into well 
separated patches after the inclusion of SNO results. The most serious consequence is 
however the considerable reduction on the significance of the SMA region. From the 
analysis of the four global rates a good part of this region becomes acceptable only at 
some marginal level. In these and next plots, the region above the line at high Am^ is 
excluded at 99% CL from the negative results of CHOOZ and Palo Verde |5^. Note 
that regions for very large Am^ ~ 10^'^ eV ^ are excluded at the 2 — 3a level without 
need of the CHOOZ result once we included the SNO data. 

The reduction of the SMA region can be explained from the tables of deviation 
(right part in Tables (0-^). Before the introduction of the CC-SNO global rate in the 
computation and minimization the residuals for the CC-SNO data corresponding to 
the best fit are very high (Table (^): 2.5a for the SMA and also for the LOW solution. 
The CC-SNO residual values for remaining regions are however smaller. When the 
CC-SNO data is incorporated into the fit the SMA and LOW regions become less 
favored. The effect is much more important for the relatively small SMA region. For 
the LOW region the initial area is large and nearby minima can be found to adjust for 
the CC-SNO data and still provide reasonable quality fits. 



3.2 Day-night SK spectra and global rates 

In this section we present the results obtained when including the SK day and night 
energy spectrum rates in addition to the total event rates of CI, Ga experiments and 
CC-SNO . 

The definition of the functions and minimization procedures used in the sta- 



tistical analysis are explained in detail in Section |2.2| . The number of experimental 
data inputs is 2 x 19 + 3 = 41. One now has four free parameters, the oscillation 
parameters Am^ and tan^ 6, the ®B neutrino flux normalization factor a and the cor- 
relation parameter 5cor- Two cases will be studied. In case A the flux normalization a 
is considered a free parameter, the number of effective d.o.f. is then 41 — 4 = 37. In 
case B, the parameter a is constrained to vary around the BPB2001 central value with 
a standard deviation a given by SSM. We now have 38 d.o.f. 

In Tables (^^ we present the best fit parameters or local minima obtained from the 
minimization of the function given in Eq. (^. As in the previous section, the values 
of 

Xmin P^^ degree of freedom {y^/n) and the goodness of fit (g.o.f.) or significance 
level of all minima are shown. In the right part of the tables, the deviations from the 
expected values or minimization residuals for a number of experimental quantities are 
listed. 

The results for the free flux minimization or case A are presented first. As a sum- 
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mary of numerical results, which appear in detail in Table (Q) and Fig. left), the 
position of the absolute minimum is located at the LMA region (Am^, tan^ ^)'^*^'^ = 
(5.4 X 10-5 eV^O.36), with ^^Jn = 0.82 (d.o.f. n=37). The significance level 
for this point is noticeably larger (g.o.f. = 82%) than the levels obtained in global 
rates only analysis. This is only in part an artifact of the statistical machinery. It 
is a satisfactory result because it basically reflects the internal consistency of the 
data: although the number of degrees of freedom have considerably increased, the 

per d.o.f is still much the same. The minimization value of the flux normaliza- 
tion is a = 1.02 while the correlation parameter is 5cor = 1-43. The next local min- 
ima are situated at the following positions: the LOW region (Am^, tan^ 6')^*^'^ = 
(1.2 X 10"^ eV^0.75), with xLn/^ = 0-85 and the two solutions at the VAC re- 
gion (Am2,tan2^)^^^ = (2.7 x 10"^° eV^ 2.70(0.32)), with xLJn = 0.85. In this 
last region the best-fit flux normalization is down to a = 0.56 while 5cor = 1-65. 
In case B, where the flux normalization is constrained to its SSM value, the results 
are very similar and are shown in Table (^ and Fig. (H, right). The position of 
the first minimum remains practically unchanged, and is located at the LMA re- 
gion (Am2,tan2^)^^^^ = (5.4 x 10'^ eV^0.38), with xLJn = 0.82 (d.o.f. n=38) 
and g.o.f.= 82%. The next local minima are: the one situated in the LOW region 
(Am^, tan^ 0)^'^^ = (7.5 x 10"^ eV^ 0.84), with x'LiJn = 0.86, and the twin solutions 
at the VAC region (Am^, tan^ ^)^^^ = (8.9 x 10"^° eV^ 1.73(0.48)), with xLJ^ = 0-9. 
Finally, in both cases A and B we have a local minimum situated at the SMA region. 
These minima are situated at Am^ ~ 7.3 x 10"^ and tan^ ^ ~ 1.3 x 10~^ with a 
considerable poorer statistical significance. 

From the Tables (|^-^|-^^ and also from landscape scans it can be seen that there 
are some non-negligible regions in parameter space, not only the best fit points. We 
are therefore justified in converting into likelihood using the expression C = e~-^^/^, 
and proceeding to study the marginalized parameter constraints. This normalized 
marginal likelihood is plotted in Figs. (0) for each of the oscillation parameters Am^ 
and tan^ 9. We present results corresponding to two cases: the global-rate analysis 
excluding CC-SNO and the full analysis including SK energy spectrum. For tan^ 6 
we observe that the likelihood function is concentrated in a region 0.2 < tan^ 6 < 1 
with a clear maximum at tan^ 6 ~ 0.5 in sharp coincidence with previous results. 
The situation for Am^ is less obvious although the region at large mass differences is 
clearly favored. The half width of these curves can be used as another estimate of the 
minimum error which can be assigned to the parameters at the present and near future 
experimental situation. 
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4 Borexino implications 



Borexino is a real-time detector for low energy (< 1 MeV ) spectroscopy. The experi- 
ment's goal is the direct measurement of the ^Be Solar neutrino flux of all flavors via 



neutrino-electron scattering in an ultra-pure scintillation liquid fS^-pS 



As before, in order to estimate the Borexino signal one convolutes the Sun-neutrino 
flux with a detector response function which includes the elastic scattering cross sec- 
tions i'e,nG~ — > l^e,^l€■~ , energy resolution and energy-dependent efficiency effects. We 
assume a Gaussian resolution function with an energy dependent width a{Eexp) = 
0.048 ^/E^p + 0.01 MeV . At this stage we suppose unit efficiency over the nominal 
window of the experiment 0.3 < E < 0.8 MeV and zero otherwise. The absolute ^Be 
flux is taken from BPB2001 as for the other experiments. After one year of data 
taking, Borexino |^ expects roughly 15 — 20 x 10^ Solar neutrino events with a nearly 
negligible statistical error astat ~ 0.8% whereas the error in discriminating signal from 
background is much larger a ~ 8% for the same period. 

In Table and Fig. (|^) we present the expectations for the normalized day- 
night signal 5''='«J/-"*9'^* at the Borexino experiment for all of the local minima solutions 
found in the previous section (SK spectrum plus global rates with constrained flux): 
^Bor _ s^~'^ (Am'^ ,9)/ So where 5*0 is the expected signal in absence of oscillations 
and day-night signals are averaged. We also present the expected day-night asymmetry 
A^'^ = 2{D — N) / {D + N). As we will see below, the asymmetry on the day-night event 
rates is a valuable tool for distinguishing among the different oscillation solutions. In 
Fig. (^) we present the expectations for the Borexino experiment as a function of the 
two-dimensional oscillation parameters. The expected signal varies mildly in ample 
regions of the parameter space. We observe that at the best flt solution, situated 
at the LMA region, the expected Borexino normalized signal is 5*^°^ = 0.63. In the 
whole LMA region (99.7% CL) the signal varies between 0.5 and 0.7. If we restrict 
ourselves to the 90% CL allowed region around the absolute minimum the signal is 
always within the range ~ 0.6 — 0.7. The asymmetry at the absolute minimum is 
negligible A^'^ < 10~^, while in the LMA region which surrounds it, it can reach the 
~ 1% level. A similar behavior appears in all the rest of the parameter space except 
in the LOW region. 

In the LOW region the signal varies between ~ 0.6 — 0.7 (90% CL). The variation 
increases to cover the range ~ 0.5 — 0.8 at 99.7 CL. Here the day-night asymmetry 
is expected to be at its maximum. From the table one observes values as high as 

ADN ^ 30%. 

In the VAC region the signal varies between ~ 0.7 — 0.8. The variation of the 
expected signal in the SMA region is much more abrupt. This prevents us from giving 



11 



accurate predictions for this case. The signal passes from a value 5''^°'' ~ 0.2 to 
5*^°^ ~ 0.8 in a very narrow region. As it is shown in Fig. (|^), the combined Borexino 
measurements of the total event rate with an error below ±7 — 10% and day-night total 
rate asymmetry with a precision comparable to that of SK will allows us to distinguish 
or at least to strongly favor the Solar neutrino solutions provided by present data. 



5 Conclusions 

We have analyzed experimental evidence from all the Solar neutrino data available at 
this moment including the latest charged current SNO results. For the best solutions 
which fit the present data, we have obtained and analyzed the expectations in Borexino. 
We considered different combinations of non-redundant data to perform two kinds of 
analyses: in the first one we only included total event rates for all the experiments, 
while in the second we included global rates for Homestake, GNO-GALLEX-SAGE 
and CC-SNO plus the day and night energy spectrum rates provided by SK. In this 
last case we studied the possibility of free or constrained variations of the Boron Solar 
neutrino flux. 

In the simplest case involving global rates only we have observed how the inclusion 
of the CC-SNO data causes the strong decrease of the statistical significance of the 
SMA solution. 

In the most comprehensive case, global rates plus spectrum, the best fit was ob- 
tained in the LMA region with parameters (Am^, tan^ q^lma ^ (^^^^ ^ j^g-s ^y2^ 0.38), 
{xl^Jn = 0.82, d.o.f. n=38). Solutions in the LOW and VAC regions are still possible 
although much less favored. The best possible solution in the SMA region gets a low 
statistical significance. 

We have analyzed in detail the expectations on the future experiment Borexino 
for all the favored neutrino oscillation solutions. The expected Borexino normalized 
signal is S^°^ = 0.63 at the best-fit LMA solution while the day-night asymmetry is 
negligible A^^ < 10^^. In the VAC region the signal is slightly higher ~ 0.7 — 0.8 
and the asymmetry is still practically negligible. In the whole LMA region (99.7% CL) 
the signal varies between 0.5 and 0.7. In the LOW region the signal is in the range 
~ 0.6 - 0.7 at 90% CL while the asymmetry A^'^ ~ 1 - 30%. We conclude that in the 
near future, after 2-3 years of data taking, the combined Borexino measurements of 
the total event rate with an error below ±5 — 10% and day-night total rate asymmetry 
with a precision comparable to that of SK will allows us to distinguish or at least 
to strongly favor the Solar neutrino solutions provided by present data. Additional 
information and precision measurement should be obtained also from Kamland and 



from future solar neutrino experiments |56 



12 



Acknowledgments 



It is a pleasure to thank R. Ferrari for many enlightening discussions and for his 
encouraging support without which this work would have not been possible. We thank 
all the Borexino group of Milano University and especially M. Giammarchi, E. Meroni 
and B. Caccianiga for providing us essential information about the experiment and its 
signal discrimination power. One of us (V.A.) would like to thank M. Pallavicini for 
interesting and useful discussions about Borexino features and potentiality. We thank 
also all those who sent us comments on early versions of this work. We acknowledge 
the financial support of the Italian MURST, the Spanish CYCIT funding agencies and 
the CERN Theoretical Division. The numerical calculations have been performed in 
the computer farm of the Milano University theoretical group (THEOS). 



13 



References 



R. Davis, Prog. Part. Nucl. Phys. 32 (1994) 13. B.T. Cleveland et al., (HOMES- 
TAKE Coll.) Nucl. Phys. (Proc. Suppl.)B 38 (1995) 47. B.T. Cleveland et al., 
(HOMESTAKE Coll.) Astrophys. J. 496 (1998) 505-526. 

Y. Fukuda et al. [Super-Kamiokande Collaboration], Phys. Rev. Lett. 82, 2430 
(1999) ||arXiv:hep-ex/9812011|| . 



Y. Fukuda et al. [Super-Kamiokande Collaboration], Phys. Rev. Lett. 82, 1810 
(1999) ||arXiv:hep-ex/9812009|| . 



J.N. Abdurashitov et al. (SAGE Coll.) Phys. Rev. Lett. 83(23) (1999)4686. 

P. Langacker. Talk given at 4th Intl. Conf. on Physics Beyond the Standard Model, 
Lake Tahoe, CA, 13-18 Dec. 1994, |hep-ph/9503327| ; Published in Trieste HEP 



Cosmol. 1992:0487-522. Ibid., Nucl. Phys. (Proc. Suppl.)B 77 (1999) 241, |hep 
ph/9811460[ Ibid.., Talk given at 1th Int. Workshop on Weak Interactions and 



neutrinos (WIN99), Cape Town, SA, 24-30 Jan 1999, |hep-ph/ 9905428 



P. Ramond, |hep-ph/980940l| , Nucl. Phys. (Proc. Suppl.)B 77 (1999) 3. 
F. Wilczek, |hep-ph/ 98095091 Nucl. Phys. (Proc. Suppl.)B 77 (1999) 511. 



S.M. Bilenky et al. Summary of the Europhysics neutrino Oscillation Workshop 
Amsterdam, The Netherlands, 7-9 Sep 1998, |hep-ph/990625Tl . S.M. Bilenky, C. 
Giunti and C.W. Kim, |hep-ph/990246^ . Int. J.Mod.Phys.A15(2000) 625. 



E. Torrente-Lujan, |arXiv:hep-ph/ 9902339 



P. Aliani, V. Antonelli, M. Picariello and E. Torrente-Lujan, arXiv:hep^ 



ph/0112101 



G. Fiorentini, F. L. Villante, B. Ricci, ||arXiv:hep-ph/0109275|] . 



Q. R. Ahmad et al. [SNO Collaboration], Phys. Rev. Lett. 87 (2001) 071301 
| ^rXlv:nucl-ex/0l060l5|| . 

[13] S. Turck-Chieze, Nucl. Phys. Proc. Suppl. 91 (2001) 73. 



E. G. Adelberger et al. Rev. Mod. Phys. 70 (1998) 1265 ||arXiv:astro-ph/9805121 
A. S. Brun, S. Turck-Chieze and P. Morel, |arXiv:astro-ph/9806272 



14 



[14] J. N. Bahcall, M. H. Pinsonneault and S. Basu, Astrophys. J. 555, 990 (2001) 
[|arXiv:astro-ph/00103l6[ . 



[15 
[16 

[17; 

[is: 



[19 



[20 

[21 
[22 
[23 
[24 

[25 
[26 

[27; 
[28 
[29 



J.N. Bahcall and M.H. Pinsonneault, Rev. Mod. Pliys. 67 (1995) 781. 

J. R. Klein [SNO Collaboration], In *Vemce 1999, Neutrino telescopes, vol. 1* 
115-125. 

A. B. McDonald [SNO Collaboration], Nucl. Phys. Proc. Suppl. 77 (1999) 43. 

J. Boger et al. [SNO Collaboration], Nucl. lustrum. Meth. A 449 (2000) 172 
||arXiv:nucl-ex/9910016| |. 



V. Barger, D. Marfatia and K. Whisnant, Phys. Lett. B 509 (2001) 19 [|arXiv:hep- 
ph/0l04l66| 



J. N. Bahcall, P. I. Krastev and A. Y. Smirnov, JHEP 0105 (2001) 015 ||arXiv:hep- 
ph/0103179|| . 



Y. Fukuda et al. [Super-Kamiokande Collaboration], Phys. Rev. Lett. 81, 1158 
(1998) [Erratum-ibid. 81, 4279 (1998)] ||arXiv:hep-ex/9805021|' . 



A. W. Poon [SNO Collaboration], |arXiv:nucl-ex/0110005 



M. V. Garzelli, C. Giunti |arXiv:hep-ph/0111254 . 

V. Barger, D. Marfatia and K. Whisnant, |arXiv:hep-ph /0 106207 



G. L. Fogli, E. Lisi, D. Montanino and A. Palazzo, Phys. Rev. D 64 (2001) 093007 
||arXiv:hep-ph/0106247 . 



P. I. Krastev and A. Y. Smirnov, |arXiv:hep-ph /0 108 1 77 . 



J. N. Bahcall, M. C. Gonzalez-Garcia and C. Pena-Garay, JHEP 0108, 014 (2001) 
[|arXiv:hep-ph/0106258|| ; J. N. Bahcall, M. C. Gonzalez-Garcia and C. Pena-Garay, 
arXiv:hep-ph/0111150 



A. Bandyopadhyay, S. Choubey, S. Goswami and K. Kar, |arXiv:hep-ph /Ol 10307 
S. Choubey, S. Goswami and D. P. Roy, |arXiv: hep-ph /0109017 . 



A. Bandyopadhyay, S. Choubey, S. Goswami and K. Kar, Phys. Lett. B 519 (2001) 
83 [|arXiv:hep-ph/0106264|| . 



[30] M. Ahmann et al. (GNO Coll.) Phys. Lett. B490 (2000) 16-26. 



15 



[31] P. Anselmann et al, GALLEX Coll., Phys. Lett. B 285 (1992) 376. W. Hampel 
et al., GALLEX Coll., Phys. Lett. B 388 (1996) 384. T.A. Kirsten, Prog. Part. 
Nucl. Phys. 40 (1998) 85-99. W. Hampel et al., (GALLEX Coll.) Phys. Lett. 
B 447 (1999) 127. M. Cribier, Nucl. Phys. (Proc. Suppl.)B 70 (1999) 284. W. 
Hampel et al, (GALLEX Coll.) Phys. Lett. B 436 (1998) 158. W. Hampel et 
al., (GALLEX Coll.) Phys. Lett. B 447 (1999) 127. 

[32] T. Hagner [BOREXINO Collaboration], Part. Nucl. Lett. 104, 90 (2001). 

S. M. Bilenky, T. Lachenmaier, W. Potzel and F. von Feilitzsch, |arXiv:hep^ 



ph/0109200i 



[34] C. Arpesella et al. [BOREXINO Collaboration], |arXiv:hep-ex/0109031 



[35] E. Meroni, Nucl. Phys. Proc. Suppl. 100, 42 (2001). 

[36] E. Torrente-Lujan, Phys. Rev. D 59 (1999) 093006. E. Torrente-Lujan, Phys. 
Rev. D 59 (1999) 073001. E. Torrente-Lujan, Phys. Lett. B 441 (1998) 305. 
V.B. Semikoz, E. Torrente-Lujan, Nucl. Phys. B 556 (1999) 353. E. Torrente- 
Lujan,Phys. Lett. B494 (2000) 255 ||hep-ph/9911"458| . 

[37] I. Mocioiu and R. Shrock, Phys. Rev. D 62 (2000) 053017 ||arXiv:hep-ph/0002149|| . 

A. Dziewonski, in The Encyclopedia of Solid Earth Geophysics, edited by D.E 
James (Van Nostrand Reinhold, New York 1989). 

[38] J. Bahcall, M. Kamionkowski, A. Sirlin, Phys. Rev. D51 (1995) 6146. 

[39] M. Passera, Phys. Rev. D 64 (2001) 113002 ||arXiv:hep-ph/0011190 |. 

[40] S. Nakamura, T. Sato, V. Gudkov and K. Kubodera, Phys. Rev. C 63 (2001) 
034617 ||arXiv:nucl-th/0009012 |. 

[41] H. Ishino, Ph. D. thesis. University of Tokio,1999. 

M. Nakahata et al. (SK Coll.), Nucl. lustrum. Methods 46 (1998) 301. 

[42] M. Nakahata et al. [Super-Kamiokande Collaboration], Nucl. lustrum. Meth. A 
421, 113 (1999) ||arXiv:hep-ex/9807027|l . 

[43] N. Sakurai, Ph.D. Thesis, Dec. 2000 Constraints of the neutrino oscillation pa- 
rameters from 1117 day observation of soloar neutrino day and night spectra in 
Super-Kamiokande. 

[44] G. L. Fogh and E. Lisi, Astropart. Phys. 3 (1995) 185. 



16 



[45] S. Goswami, D. Majumdar and A. Raychaudhuri, Phys. Rev. D 63 (2001) 013003 
[|arXiv:hep-ph/ 00031631 . 

[46] S. Fukuda et al. [Super-Kamiokande Collaboration], Phys. Rev. Lett. 86, 5656 
(2001) [|arXiv:hep-ex/0103033|] . 

[47] S. Fukuda et al. [SKamiokande Collaboration], Phys. Rev. Lett. 86, 5651 (2001) 
||arXiv:hep-ex/0103032i . 

[48] A.I. Abazov et al. (SAGE Coll.), Phys. Rev. Lett. 67 (1991) 3332. D.N. Ab- 
durashitov et al. (SAGE Coll.), Phys. Rev. Lett. 77 (1996) 4708. J.N. Ab- 
durashitov et al., (SAGE Coll.), Phys. Rev. C60 (1999) 055801; |astro-ph/990TT3T 



J.N. Abdurashitov et al., (SAGE Coll.), Phys. Rev. Lett. 83 (1999) 4686; |astro^ 
ph/9907113| 



[49] D. E. Groom et al. [Particle Data Group Collaboration], Eur. Phys. J. C 15 (2000) 
1. 

[50] M. ApoUonio et al. (CHOOZ coll.), |hep-ex/ 990703^ , Phys. Lett. B 466 (1999) 
415. M. Apollonio et al, Phys. Lett. B 420 (1998) 397. 

F. Boehm et a;.,Phys. Rev. D62 (2000) 072002 ||[iep-ex/ 000302^ . 



[51] J. Burguet-Castell and O. Mena, [axXiv:hep-ph/010"8T0£ . 

[52] K. Lande (For the Homestake Coll.) Nucl. Phys. B(Proc. Suppl.)77(1999)13-19. 

[53] J.N. Bahcall, P.I. Krastev and E. Lisi, Phys. Rev. C 55, 494 (1997); 

B. Faid, G.L. Fogli, E. Lisi and D. Montanino, Astropart. Phys. 10, 93 (1999). 



[54] V. Berezinsky, |arXiv: hep-ph /0108166 



[55] M. Giammarchi (Borexino Coll.), private communication 

A. Strumia, F. Vissani, JHEP 0111 (2001) 048 ||arXiv:hep-ph/0109172|| ; A. De 
Gouvea, C. Peiia-Garay, Phys. Rev. D 64 (2001) 113011 [ |arXiv: hep-ph/ 0107186 ] 



17 



Experiment [Ref.] 


SsSM 


Soata/ SsSM (±lcr) 


SK (1258d) H 
CC-SNO (240d, 8-20 MeV) |12| 
CI [0 
SAGE [|, 48| 
GNO-GALLEX 31] 


2.32 ±0.03 ±0.08 
2.39 ±0.34 ±0.15 
2.56 ±0.16 ±0.16 
67.2 ±7.0 ±3.2 
74.1 ±6.7 ±3.5 


0.451 ±0.011 
0.347 ±0.029 
0.332 ±0.056 
0.521 ±0.067 
0.600 ±0.067 



Table 1: Summary of data used in this work. The expected signal (Sssm) and observed 
ratios Soata/ Sssm with respect to the BPB2001 model are reported. The SK and CC-SNO 
rates are in 10^ cm^^ s^^ units. The CI, SAGE and GNO-GALLEX measurements are 
in SNU units. In this work we use the combined results of SAGE and GNO-GALLEX: 
Seal Sssm {Ca = SAGE+GALLEX+GNO)=0.579 ± 0.050. The total flux is taken 
from the BPB2001 model ||: (t>u{^B ) = 5.05(ll[j;iJ]) x 10^ cm'^ s"^ . 



Solution 










g.o.f. 


SK 


CI 


Ga 


SNO 


ADN 


SMA 


7.8 xlQ- 


-6 


1.10 xlO-^* 


1.1 


29 


0.5 


0.8 


0.2 


2.5 


1.3 


VAC 


8.4 xlO- 


-11 


0.27 


1.3 


26 


1.9 


1.4 


0.8 


0.2 


1.3 








3.37 


1.6 


21 


1.7 


1.4 


0.9 


0.3 


1.3 


LMA 


1.2 xlQ- 


-5 


0.34 


1.9 


17 


0.0 


1.3 


0.2 


0.7 


1.1 


LOW 


1.1 xlO- 


-7 


0.75 


3.6 


5.7 


2.0 


1.9 


1.0 


2.5 


0.3 



Table 2: (Left) Best fit oscillation parameters: Am^(eV^) and ia,v?6, from the analysis of 
global rates: CI, Ga and SK. The SNO measurement is not included (d.o.f.= 3 — 2). The 
significance level g.o.f. is in percentage. The sampling error is ~ 1% for Am^ and tan^(^). 
(Right) Minimization residuals in la units for diverse experimentally measured quantities: 
the global rates for the four experiments and SK global day-night asymmetry {A^^). Only 
the first three quantities appear in the computation. 



Solution 






tan^e 




g.o.f. 


SK 


CI 


Ga 


SNO 


ADN 


LMA 


3.1 xlO^ 


-5 


0.40 


1.6 


21 


1.4 


0.9 


1.5 


1.3 


1.6 


SMA 


7.3 xlO- 


-6 


1.40 xlO-^* 


2.4 


8.8 


1.0 


1.5 


0.4 


0.9 


1.3 


VAC 


9.7 xlO- 


11 


0.26 


3.1 


4.6 


3.8 


2.3 


1.1 


1.5 


1.3 








4.00 


3.2 


4.1 


4.1 


2.5 


1.3 


1.9 


1.3 


LOW 


3.2 xlO- 


-7 


0.80 


3.3 


3.5 


2.7 


1.2 


3.1 


1.7 


0.9 



Table 3: Best fit oscillation parameters, Am^(eV^) and tan^ ^, and minimization residu- 
als (see explanation in Table ^). The analysis now includes the global rates for the four 
experiments: CI, Ga, SK and CC-SNO (d.o.f. =4-2). 
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Solution 






tan^e 




g.o.f. 


SK 


CI 


Ga 


SNO 


ADN 


LMA 


5.4 xlO- 


-5 


0.36 


0.8 


82 


1.8 


0.4 


0.8 


2.1 


0.0 


LOW 


1.2 xlO" 


-7 


0.75 


0.8 


80 


2.0 


2.0 


0.9 


2.5 


0.2 


VAC 


2.7 xlO- 


-10 


2.70 


0.8 


80 


2.2 


1.8 


0.5 


0.7 


1.2 








0.32 


0.9 


80 


2.3 


2.0 


0.6 


0.8 


1.2 


SMA 


7.3 xlO- 


-6 


1.30 xlO-^' 


1.1 


40 


0.1 


1.0 


0.3 


2.1 


1.3 



Table 4: Best fit oscillation parameters, Ar?T,^(eV^) and tan^0, and minimization residuals 
(see explanation in Table |2|). The analysis now includes the global rates for three experiments 
Cl,Ga and CC-SNO , and SK day and night energy spectra. There are four parameters: 
Am^, tan^ 9, a and 5cor- We let the flux normalization a vary freely, so the d.o.f is = 41 — 4 
(see the text). The sampling error due to the finite-grid is ^ 1% for Am^^„ and tan^(0). 



Solution 


Am^ 




tan^ e 




g.o.f. 


SK 


CI 


Ga 


SNO 


ADN 


LMA 


5.4 xlO- 


-5 


0.38 


0.8 


82 


1.6 


0.5 


0.9 


1.8 


0.0 


LOW 


7.5 xlO- 


-8 


0.84 


0.9 


72 


2.7 


2.3 


0.8 


2.9 


0.7 


VAC 


8.9 xlO- 


-10 


1.73 


0.9 


73 


2.4 


1.9 


0.5 


0.7 


1.3 








0.48 


0.9 


63 


2.4 


2.2 


1.4 


1.6 


1.3 


SMA 


7.3 xlO- 


-6 


1.30 xlO-^ 


1.1 


40 


0.1 


1.0 


0.3 


2.1 


1.3 



Table 5: Same as Table Q except that here the flux normalization a is constrained to vary 
with the SSM standard deviation and the d.o.f is = 41 — 3 (see the text). 



Solution 


Am2(eV2) 


tan2(6l) 




(%) 


LMA 


5.4 xfO-5 


0.38 


0.63 


-0.008 


LOW 


7.5 xlO-s 


0.84 


0.61 


-14.3 


VAC 


8.9 xlO"i° 


1.73 


0.64 


+0.02 






0.48 


0.65 


-0.07 


SMA 


7.3 xlO-s 


1.30 xlO-3 


0.24 


-0.02 



Table 6: Borexino experiment expectations for the normalized total signal {S^"^) and day- 
night asymmetry (^4^^) for all of the best-fit solutions appearing in Table (|5|). 
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Global: Pre-SNO 



Global: Post-SNO 
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1-4 




10' 



,0 



tan2 e 



tan2 e 



Figure 1: Best fit solutions for the global rate analysis. Global rates from SK, CI and Ga 
experiments are included. The right figure contains also the CC-SNO one. The black dots 
are the best fit points (the absolute minimum is located at the SMA (left) or LMA (right) 
regions). The colored areas are the allowed regions at 90, 95, 99 and 99.7% CL relative to 
the absolute minimum. The region above the solid line is excluded by CHOOZ results at 



99% CL g. 
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tan2 e tan2 6 



Figure 2: Global solutions for the analysis of C/, Ga and CC-SNO experiments plus the 
day-night SK spectrum rates. The black dots correspond to best fit points. The colored areas 
are the allowed regions at 90, 95, 99 and 99.7% CL relative to the absolute minimum. Left 
(case A of the text): the minimization with respect 5cor is performed as before. The absolute 
normalization factor a is allowed to vary freely. Right (case B of the Text): the absolute 
normalization factor a and correlation parameter Jcor are varied subject to constraints. 




Am2 /eV2 



Figure 3: Marginalized likelihood distributions for each of the oscillation parameters Am^ 
(right), ia-n^e (left) appearmg m the fit- The curves are in arbitrary units with normal- 
ization to the maximum height. The continuous lines shows the distributions corresponding 
to Fig. dll, right). Dashed lines correspond to the case represented in Fig. (||, left). 
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tan2 e 

Figure 4: The signal at the Borexino experiment as a function of the oscillation parameters. 
The signal is normalized to the no-oscillation case. Contours are drawn at 5/5*0 = 0.5, 0.6, 0.7 
with full lines (respective from inside to outside). Superimposed: the allowed regions from 
the global rate analysis including CC-SNO (see Fig. (^), right) and the la allowed regions 
(here enclosed by the dashed lines) from the ^Be sensitive CI experiment alone. 
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0.2 0.4 0.6 0.8 1 



Figure 5: Day night asymmetry | A \ versus normalized signal at Borexino corresponding 
to the 90, 95, 99, 99.7 % CL allowed regions from the full analysis including global rates and 
SK spectra. The upper (lower) zones correspond to the LOW and VAC (LMA) regions in 
Fig. (0, left). The stars are the expectations of the best fit solutions in each of the regions. 
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